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Abstract 



The motion of neutral particles with magnetic moments in an inhomoge- 
neous magnetic field is described in a semi-classical framework. The concept 
of Coherent Internal States is used in the formulation of the semiclassical 
approximation from the full quantum mechanical expression. The classical 
trajectories are defined only for certain spin states, that satisfy the conditions 
for being Coherent Internal States. The reliability of Stern-Gerlach experi- 
ments to measure spin projections is assessed in this framework. 
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I. INTRODUCTION 



The Stern-Gerlach experiment consists in taking a beam of particles that have a neutral 
electric charge, but a finite magnetic moment, and making them to go through an inhomoge- 
neous magnetic field. The observed result is that the particles deflect differently depending 
on the spin projection along the magnetic field. So, by measuring the deflection, one can 
infer the value of the spin projection of the particles along the direction of the magnetic 
field. The Stern-Gerlach experiment is the archetype of the measurement of a quantum me- 
chanical property. Thus, it is always discussed even in the most basic textbooks of quantum 
mechanics []IJ]. The explanation that it is usually done for the Stern-Gerlach experiment is 
of a semiclassical nature. The motion of the particles is approximated by classical trajecto- 
ries. In the following, we will present a brief account of the description of the Stern-Gerlach 
experiment, as it is usually done in textbooks. 

One starts with a beam of particles, in a certain spin state, moving initially along a 
straight line. For the following discussion, we will take the y-axis along the direction of the 
motion of the particles. The center of the beam corresponds to the coordinates x = 0, z = 0. 
The beam enters in a Stern-Gerlach magnet. Usually, the Stern-Gerlach magnets produce 
magnetic fields that are independent of y and that do not have components in the y-direction 
(neglecting border effects). The z-axis is chosen along the direction of the magnetic field at 
the center of the beam. So, 

B(f) = B z (x,z)u z + B x (x,z)u x (1) 

where B x (0, 0) = 0. The geometry of the Stern-Gerlach magnet is such that B z (x, z) varies 
with z, but is mostly independent of x. Thus, a force appears in the z-direction which is 
proportional to dB z /dzfi z . However, as the divergence of the magnetic field has to vanish, 
then dB z /dz = —dB x /dx, and so if the B z varies with z, then B x must vary with x. Thus, 
if we only retain terms linear in x and z, we get 

B(r) = {B -B x z)u z + xB x u x (2) 

The magnetic field interacts with the magnetic moment of the particle fi, giving rise to 
an interaction energy given by 

E B {r) = -B(r) ■ ft (3) 

This interaction energy depends on the distance. Thus, a force is generated which is given 
by 

F = -VE B (r) = VB z (x, z)fi z + VB x (x, z)fjL x (4) 

Thus, a force in the x-direction also appears which would be proportional to dB x j dx[i x . 

^From the expression ([|), the force acting on a particle with a certain spin state can be 
calculated. If we consider the scattering of particles which have a given spin projection m 
along the z-axis, they will suffer a force that will be calculated as the expectation value in 
the z-direction which is proportional to the spin projection. The second term in (f4|) will not 
contribute because the matrix element < m\fi x \m > vanishes. This result is well known, 
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and it is in agreement with experiment. However, we can consider the case of the scattering 
of particles with spin projection ml along the x-axis. If we calculate the expectation value 
of F in (|]), we obtain a force in the x-direction, while the force in the z-direction vanishes. 
This implies a deflection of the trajectory in the x-direction, that does not happen in the 
experiment. As any student of quantum mechanics should know, the trajectory is split in 
21+1 trajectories, each one of which is deflected in the z-direction by a different amount. 

In most of the basic textbooks of quantum-mechanics, the term proportional to fi x in 
is simply ignored . Other more recent books, such as || and || , only consider the spatial 
variation of the component B z of the magnetic field, but they neglect the variation of B x . 
The classic book of quantum mechanics by Messiah argues that fi z is basically constant, 
while fi x oscillates around zero. So, in the average of the force over many oscillations this 
term would cancel. This argument, although plausible, is hardly a firm ground on which to 
describe the general motion of particles in inhomogeneous magnetic fields. The opinion of 
the authors is that the Stern-Gerlach experiment still requires a satisfactory description in 
semiclassical terms. 

In this paper we propose a semiclassical explanation to the Stern-Gerlach experiment. 
It relies on the concept of Coherent Internal States (CIS), that we introduced in a previous 
paper |5|]. The conclusion of that paper was that when one analyzes the scattering of a 
particle with internal degrees of freedom, such as the spin, a single trajectory is a meaningful 
approximation for the quantum mechanical scattering wave function only for a certain set 
of internal states, which we called Coherent Internal States (CIS). Thus, if the scattered 
particle has an internal state that coincides initially with one of the CIS, then its scattering 
wavefunction can be approximated by a single trajectory. If not, the internal state should 
be expanded in terms of the CIS, and then the scattering wavefunction can be approximated 
by a combination of classical trajectories, one for each CIS. 

In the case of Stern-Gerlach experiments, we demonstrate in section 2 that the CIS 
correspond to states with definite projection along the direction of the magnetic field. This 
direction, called z', may vary depending on the position of the particle, because the magnetic 
field is not homogeneous. So, it does not coincide with the laboratory fixed z-axis defining 
the direction of the magnetic field at the centre of the beam. In section 3 we use the path 
integral formalism to describe the trajectories of the CIS, and we find that they deflect 
on the z' direction. Considering that the average of z' corresponds to z, that explains the 
observed fact that the trajectories split in the z-direction, and not in the x-direction, as 
(|) could suggest. In section 4 we describe the scattering of a beam of particles with a 
finite size in a Stern-Gerlach magnet. As the axes z' and z do not coincide, the deflection 
of the trajectory is not always consistent with the spin projection along the z axis. So, a 
Stern-Gerlach magnet, understood as a measurement apparatus to find the spin projection 
along the z axis, has a certain probability of giving a wrong result, which can be evaluated 
in our formalism. Section 5 is for the summary and conclusions. 

II. COHERENT INTERNAL STATES FOR A PARTICLE MOVING IN AN 
INHOMOGENEOUS MAGNETIC FIELD 

It was shown in || that the notion of a classical trajectory is a useful approach for the 
quantum mechanical wave function only for certain selected states, that were called Coherent 
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Internal States (CIS). These states are the eigenstates of the cross section matrix, which are 
orthogonal and form a basis of the space of the initial internal states. To find the CIS, as it 
is was shown in ||, the following iterative procedure should be followed: 

• Solve the classical scattering problem for the uncoupled hamiltonian H and obtain 
the evolution operator along the classical trajectory. 

• Consider small desviations from the classical trajectory. Evalute the operator A, 
defined in |9| , which describes the dependence of the cross section matrix on the initial 
state. 

• Obtain the CIS \n > diagonalizing the cross section matrix, which is equivalent to 
diagonalize the operator A. 

• Evaluate the classical trajectories r n (t) for each CIS, the evolution operators U n (t,t ), 
and the final states |n >. If the final states are orthogonal, then the calculated cross 
section matrix will be diagonal, and the self-consistency would have been achieved. If 
not, the CIS should be recalculated as the eigenstates of the cross section matrix and 
the procedure should be followed until self-consistency is achieved. 

Let us consider a neutral particle of mass M that moves in the y-direction with velocity 
v y and which has an initial position characterized by the coordinates (xq,Zo). Note that 
in strict quantum mechanical terms, we can take wave packets sufficiently localized around 
(xq,Zq), which would have momentum dispersions much smaller than Mv y . This particle, 
that has a spin / and a magnetic moment /J = hqI, enters in a magnetic field. The time 
evolution of this particle moving in a magnetic field is given by the evolution operator 
U(tf,to) = exp{—i/h f t * Hdt}, where the Hamiltonian H of the system is written as 

2 

H = ^tj - At- B (5) 

The first term is the kinetic energy of the particle and the second term is the potential 
interaction of the particle with the magnetic field. The magnetic field has x, z components 
so that, 

pi-B = fi [(B - B x z)I z + xB 1 I x ] (6) 

We will define new axes (x', z') so that the magnetic field B(x , z ) is directed along z'. 
The y axis is unaffected. The angle ft(x , z ) that generates the rotation is given by 

tan(/3(x , zo)) = x B 1 / (B - zqB x ) (7) 

The angular momentum operators in the new coordinate system are 

I' z — Iz cos (3 + I x sin p 
I'x = -Iz sin /3 + I x cos (3 (8) 



The interaction term with the magnetic field is given by, 
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ft - B = HqB'qI's + hqB x [-(z - z )I z + (x- x )I x 



(9) 



with 



B' = tJ(B - B lZQ f + (^xo) 2 
The Hamiltonian (||) can be written as: 



H = H + V 

,2 



Ht 



o 



P 
2M 



fj, B' r 



V = (MoB^-iz - z )I z + (x- x )I x 



(10) 

(11) 
(12) 
(13) 



For the case that we are considering, the classical trayectories for H are just straight lines, 
given by the expressions: 



r( ) = ro+ M /_( / " } ^ 



M 



(14) 



We consider the trajectory of a particle that is initially in the position and 
moves in the y direction. The evolution operator for H is: 



U {t,t ) = exp{-[(t - t )n B' o r z )}. 



(15) 



In a basis {|m >} of eigenstates of I' z , the matrix elements of the evolution operator associ- 
ated with the classical trajectory are diagonal: 



< m f \T°\mi >=< mf\U (tf,t )\mi >= S mf ^ exp{-iu (t f - t )rni} 



(16) 



with u = fi B' /h. 

When the effect of small deviations from the classical trajectory in the path integral 
formalism are considered, the expression for the scattering amplitude becomes 



< mf\T\mi > = < rrif\T \rrii > + < mf\8T\mi >= 

< rrif\T \rrii > + < rrif\T x \mi > + < m/|T'|mj > 



(17) 



The correction terms are given by 

1 ftf 



< mf\T l \mi > 

< mf\T'\mi > 



2 J t 



dt < rrif\Uo(tf, t) 







dfn I - \ 

u / po ^ 



( dV(t) \ 



2h Jto 



to 



V dr J 



dt I dt' < m f \U (t f ,t) 
U (t,t'] 



V d Po J 



U (t,t )\rn i > (18) 



PO 



I'll 



f dVjf) 
V dp 



'■(I. 



Ua(t',to)\mi > 



(19) 



The derivatives w.r.t. po can be expressed in terms of derivatives w.r.t. f using eq. ([L4] 
The term < m/lT^mj > vanishes because V 2 B = and so V 2 V(t) = 0. The term T' m i 
can be calculated in a straightforward way resulting: 
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<m f \T'\m i >= OoSi) 2 ^ jf ' dt £ dt'(t-t') <m f \[U {t f ,t)I x U {t,t')I x U {t',t ) 

+U (t f , t)I z U (t, t')I z U (t', t )] \rrn > (20) 

This expression can be written in terms of the rotated angular momentum operators I' x , I' z : 

< m f \T'\m t >= (/io^) 2 ^ ttfMtt ~ < m f \ [U (t f ,t)I' x U (t,tf)I' x U (lf,t ) 

+U (t f ,t)I f z U (t,t')I' z U (t',t )}\m i > (21) 

Then, the first-order correction A m / m defined as; 

< m'^rrii >= £[< 77i/|T°K >*< m f \5T\mi > + < m/lT^m* >< ™/l<^K >*] (22) 
is given by: 

<m' i \A\m i >= (/io^i) 2 -^- f dt\t-t') < m\\ \u (t , t)I' x U (t, f)I' x U (f,t )+ 

ZTbM Jto Jt 1 

U (t ,t)r z U (t,t')r z U (t',t )} \rrii > +h.c. (23) 

The non-diagonal matrix elements of (23) vanish. The diagonal matrix elements are given 
by: 

< m-|A|mi >= Sm' m. 2r "° ^3 ° [2 - 2cos(w (t/ - *o)) - (u (t f -t )smu (t f -t ))] (24) 

ivjf /I Ct^o 

This result shows that the states \rrii >, with a definite spin projection along the z' 
axis, are the eigenstates of A, and thus they are our initial choice for CIS. As we will see 
in the next section, these states are not modified as the particle moves along the classical 
trajectory, and so they are the Coherent Internal States of our problem. 

The original states \m; L >, which have a definite spin projection m along the laboratory 
fixed z-axis, are not CIS. To describe the evolution of these states, they should be expanded 
in terms of the CIS, by means of the expression: 

\m;L>=J2d I m , tm ((3(x ,z ))\m' > (25) 



III. SEMICLASSICAL DESCRIPTION OF THE SCATTERING FOR COHERENT 

INTERNAL STATES 

We will now evaluate the classical trajectory for each CIS, by making the stationary 
phase approximation on the matrix elements of the exact propagator on the calculated CIS. 

The exact propagator of the system between the CIS is given as a path integral extended 
to all possible trajectories by: 

< m(f)\K\m >= J D[r(t)} exp{^ e// (f(t))} (26) 
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where 



S eff = f -Po + S (r(t)) ~ihln< m{f)\U B {r{t);t f ,t )\m > (27) 
is the effective action, So(r(T)) is the action corresponding to the kinetic energy, and 

U B {r{t)] t f , t ) = exp{^ jf*' dtfL ■ B(f(t))} (28) 

is the evolution operator of a particle in magnetic field along the trajectory f(t). In this 
expression, \m > is one of the CIS, and \m(f) > is the final state, defined by 

\m(f) >= U B (f m (t);t f ,t )\m > (29) 

The classical trajectory f m (t) is obtained imposing the stationary phase condition in ([26]), 



' >fl f(t)=f m (t) = (30) 



5f(t) 

This leads in a strightforwad way to: 

dr m (t) po_ 

dt lt=to ~ M 

d 2 f (t) 1 B 

= ~m< rn(t)\^(f2-B(m?=? mi t)\rn(t) > (31) 

where the state at the instant t is given by \m{t) >= UB(r m {t)]t,t )\m >. These equations 
describe the motion of a classical particle with a magnetic moment /7(t) =< m(t)\fl\m(t) > 
moving in the inhomogeneous magnetic field B{f). However, it should be stressed that this 
interpretation is only meaningful for internal states \m(t) > that evolve from a CIS state 
|m > at t — to. 

We observe that in the y direction we have a constant motion given by yj = yo + (tf—to)jj 
and in the x, z directions the motion is accelerated and the force is proportional to B\ 
because: 

d 

< m(t)\ — (fl- B(r))\m(t) >= -/i 5i(< m(t)\I z \m(t) > u z - < m{t)\I x \m{t) > u x ) (32) 
or 

Known the force, the specific nature of the classical solution in the xz directions depend 
enterily of the evolution operator Ub and the \m > state. Some considerations must be 
done to solve the classical equations in the x, z directions. It should be noticed that the 
trajectory of the particle, that initially is on (xo, z ), is always along the line that joins this 
point with the point with coordinates (x = 0, z = Bq/Bi). Along this line, the direction 
of the magnetic field is fixed, although its magnitude will vary. The state \m(t) > is given 
by the initial CIS \m(t) > times a phase factor. The angle (3 between the axes z and z' is 
constant. Thus, we have 

< m(t)\I z \m(t) >= cos(f3(x ,z )) ; < m(t)\I x \m(t) >= sm(f3(x , z )) (33) 

So, the trajectory is uniformly accelerated along the x and y directions, so that, when 
they leave the magnetic field at a time tf, the coordinates have changed to 
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UqBiTTI . c, , n , , x , , 

Zf = Z ° 2M~^ f ~ to ) C0S (P( X °> z ^> ( 34 ) 

x f = x + 2M (tf-to) sm((3(x , z )) (35) 

When the particle leaves the magnetic field, the force vanishes. However, as it has acquired a 
certain velocity, the values of these coordinates at a later time td, in which they are detected, 
is given by 

z d = z - bm cos((3(x , z )) (36) 
x d = x + bmsm(/3(x , z )) (37) 

where b = /j,oB 1 (tf — to)(2td — tf — t )/(2M), measures the spacial separation between the 
different magnetic substates. Note that, by making t d sufficiently large, we can obtain 
reasonable values of b, even if the gradient of the magnetic field B\ is small, and the length 
of the magnet, which is L = (tf —t )v y , is small. 



IV. RELIABILITY OF STERN-GERLACH EXPERIMENTS TO MEASURE SPIN 

PROJECTIONS 



We will consider a beam of particles with a finite extension. Thus, the center of the 
beam will be placed at x = 0, z = 0, but it will have a probability density P(x, z) of being in 
the position (x,z). We will assume, to make the calculations simpler, that this probability 
density is gaussian: 

P(x,z) = — 2 eM-^ J ^) (38) 

7TGT GT 

so, a is a measurement of the size of the beam. Besides, we will assume that all the velocity 
of all the particles goes essentially along the y direction, with a velocity v y . That means 
that the x and z components of the velocity must be very small so that v z (td — t ) <^ b. The 
uncertainty principle implies that Mv z > h/a, and so one gets, for (td — t ) ^> (tf — 1 ), the 
condition 

B clL 

h < B ia (t f -t ) = — — . (39) 

Vy 

This condition implies that the gradient of the magnetic field cannot be arbitrarily small. 
However, we will see that it cannot be too large either. 

Let us consider the particle which is initially in the position (x,z). The magnetic field 
that it sees is given by 

B(x, z) = (B - zB x )u z + xB x u x (40) 

Thus, the direction of the magnetic field does not go along the z-direction. The angle (3(x, z) 
between B(x, z) and u z is given by 

tan f3(x, z) = xB 1 /(B - zB x ) (41) 
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Note that the states \m; L > with spin projection m along the laboratory fixed z-axis are 
no longer the CIS. For a given value of (x, z), the CIS are states \m' > which have a given 
spin projection w! along the axis z' which is paralell to B(x,z). Then, the state \m;L > 
will not give rise to a unique trajectory, but to 21 + 1 trajectories which are characterized 
by the states \m' >. The probability that the state \m; L > follows the trajectory of \m' > 
is simply 

p(m,m') = | <m;L\m' > \ 2 = |^ m ,(/3(x, z))\ 2 . (42) 

We have considered the case of a spin-1 particle. It has three possible spin projections 
along the laboratory fixed z-axis, corresponding to the states \m; L >, for m = — 1,0, +1. 
When a beam of particles with these spin states go through a Stern-Gerlach magnet, they 
will be deflected according to the value of the spin projection m' along the z'-axis. The 
probability that the deflection of the particles is determined by the spin projection along the 
laboratory fixed z-axis, is given by the expectation value of p(m, m) averaged over the beam 
distribution, is shown in figure 1. The values < p(l, 1) > and < p(0, 0) > are represented as 
a function of the dimensionless parameter aBi/B . We see that only for values aBi/B 1, 
the values of < p(l, 1) > and < p(0,0) > tend to one. This indicates that the deflection 
of the particles will be determined by the spin projection in the "laboratory" fixed axis z 
only if aBi/B Q <^ 1. Thus, this is the condition for a Stern-Gerlach experiment to be a 
measurement of the spin projection. 
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FIG. 1. Probability that a particle with spin 1 and projection m along the laboratory fixed 
z-axis is deflected according its spin projection, as a function of the parameter aBi/Bo, for m = 
and m = 1. The case m = — 1 coincides with m = 1 

If we will consider the case in which aB\ <C B , then, the angle f3(x, z) is small, and the 
expression above can be expressed as 



The Stern-Gerlach experiment is used to measure spin projections along the z-axis. Thus, 
if we observe a deflection corresponding to the state to', we would conclude that the spin 
projection along the laboratory fixed z-axis was to', instead of to. Thus, the probability 
that the measurement gives the incorrect result is 




(43) 



(44) 



(45) 




(46) 



If we average this probability over the beam probability density, we obtain 




(47) 
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Thus, we see that the crucial condition for the Stern-Gerlach experiment to be useful in 
order to measure spin projections is that the magnetic field B should be much larger than 
the product of the gradient of this field multiplied by the size of the beam. In other words, 
the relative change of the magnetic field within the finite extension of the beam sould be 
very small. Thus, we can write 

B > B 1 a. (48) 
Putting together the conditions in eqs.(^,^), we have 

MlflA » |l » i. (49) 
n B\a 

However, the first term is just the precession angle of the magnetic moment operator about 
the z-axis. This angle has to be very large, compared to 1, as a neccesary condition for eq. 



(f49|) to be valid. In this sense, our results are in agreement with the argument of Messiah 
0], which points to the fact that the x-component of the magnetic moment oscillates around 
zero, and then it can be ignored. However, we find that this argument is not sufficient. The 
gradient of the magnetic field has to be such that is much smaller than the precession 
angle, and much larger than one. 

We have performed a simulation of a beam of 1000 particles, with spin 1 and projections 
along the laboratory fixed z-axis m = —1, 0, +1, and a probability density of having initial 
(x, z) values given by equation (^). The particles go through an inhomogeneous magnetic 
field, and are deflected. The values of the magnetic field, its gradient, the length of the 
magnet, and the distance of the detectors are taken so that the parameter b, wich determines 
the amount of the deflection, is given as b = 4a, in terms of the initial size of the beam. The 
relation of the magnetic field and its gradient is given by Bia/B Q = 0.25. In figure 2 we 
present the results of the simulation, presenting the final values of the x and z coordinates, 
in units of a, for different values of m. We find that, in general, most of the particles suffer 
deflections along the z axis consistent with their spin projection, but there are a few cases 
(corresponding to about 3% for m = 0, and 1.5% for m = ±1), in which this is not the 
case. This is what we expect from figure 1. The other effect that one sees in figure 2 is a 
focusing effect for the particles with spin projection m = — 1. That is related to the fact 
that the deflection occurs along lines that cross in the point with (x = 0, z = B /Bi). The 
states with m' — — 1 tend to get close to this point, and so they focus, while the states with 
m' = 1 tend to separate from it, and so they de- focus. 
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FIG. 2. Results of a simulation of 1000 particles with spin 1 and projections m = — 1,0, +1 
which go through an inhomogeneous magnetic field described by b = 4a and Bicl/Bq = 0.25. The 
final position of the particles in the (x, z) plane is represented in units of the initial beam size a. 

In this work we have focussed in the application of Stern-Gerlach magnets as a mea- 
surement apparatus to determine the spin projection of individual atoms, assuming that the 
spin and the magnetic moment is previously known. We find that these experiments are not 
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completely reliable, because the deflection is not uniquely determined by the spin projection. 
However, Stern-Gerlach experiments may also be used to measure the magnetic moment. In 
this case, one should measure the separation of the piles of particles coming from an initially 
unpolarized beam. Our calculations could be useful for this purpose, because they not only 
give the separation, but also give the shape of the different piles. 

V. SUMMARY AND CONCLUSIONS 

We describe the motion of a particle with spin in an inhomogeneous magnetic field, 
in a semiclassical approach. We make use of the fact that the classical trajectories are 
only a meaningful approach to the quantum mechanical scattering wavefunction for certain 
states of the internal variables, that are called Coherent Internal States. The Coherent 
Internal States are obtained initially as the eigenstates of the cross section matrix. Each 
one of these states has a trajectory that describes the time dependence of the coordinate, 
and an evolution operator that describes the time dependence of the internal state. The 
trajectory and the evolution operator are related self-consistently, because the classical force 
that defines the trajectory is related to expectation value of the coupling potential on the 
internal state, while the evolution operator is related to the coupling potential evaluated 
along the trajectory. 

We have considered the case of a particle with a given magnetic moment, moving initially 
in the y-direction within a magnetic field that depends linearly on the coordinates x and z. 
The Coherent Internal States correspond initially to definite projections of the spin along 
the direction of the magnetic field evaluated at the initial point of the trajectory. The 
trajectory corresponding to a given spin proyection will be deflected due to the gradient 
in the magnetic field. The magnitude of the magnetic field observed by each particle may 
change as a result of the deflection of the trajectory. However, the direction of the magnetic 
field remains constant. Thus, the Coherent Internal State does not change as the particle 
moves along its trajectory. 

We have considered an ensambe of particles, all moving in the y-direction, with initial 
x and z coordinates following a gaussian probability distribution. This illustrates a realistic 
situation for a beam of particles of finite size entering a Stern-Gerlach magnet. We have 
evaluated the trajectories followed by these particles, considering explicitly that the Coherent 
Internal States are different for different particles of the ensamble, because they depend on 
the initial values of x and z. We find that not all the particles having a given spin projection 
along the z-axis suffer the same deflection. This indicates that even an idealized Stern- 
Gerlach experiment has a finite probability of giving the wrong result as a measurement 
apparatus of the spin projection. The probability of error depends on the relative dispersion 
of the values of magnetic field within the beam size. 

Acknowledgements: This work has been partially supported by the Spanish CICyT, 
project PB98-1111 



13 



REFERENCES 



[1] R. Eisberg and R. Resnick, Quantum Physics, Wiley, 1974. 

[2] J.M. Levy-Leblond and F. Balibar, Quantics, North Holland, 1990. 

[3] E. Merzbacher, Quantum Mechanics, Wiley, 1998. 

[4] A. Messiah, Mechanique Quantique, Dunod, 1965. 

[5] S. Cruz-Barrios and J. Gomez-Camacho, Nucl. Phys. A636 (1998) 70-84 
[6] D.M. Brink, Semi- Classical Methods for Nucleus-Nucleus Scttering, Cambridge Univer- 
sity Press, 1985. 
[7] P. Pechukas, Phys. Rev. 181 (1969) 174. 

[8] F.D. Dos Aidos, C.V. Sukumar and D.M. Brink, Nucl. Phys. A448 (1986) 333-364. 

[9] C.V. Sukumar and D.M. Brink, Nucl. Phys. A560 (1993) 863-878. 
[10] C.V. Sukumar and D.M. Brink, Nucl. Phys. A587 (1995) 413-420. 
[11] K. Alder and A. Winther, Electromagnetic Excitation, North Holland, Amsterdam, 
1975. 

[12] G. R. Satchler, "Direct Nuclear Reactions", Oxford University Press 1983, p.343. 



14 



